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§1. Introduction 

The origin of the nuclear force is one of the major unsolved problems in particle 
and nuclear physics even after the establishment of the quantum chromodynamics 
(QCD). Recently, three of the present authors proposed a new approach to extract 
the NN interactions below inelastic threshold in lattice QCDP®'® Through the 
Nambu-Bethe-Salpeter (NBS) wave function, the energy-independent but non-local 
potential U(r, r') is so defined that the NBS wave function obeys the Schrodinger 
type equation in finite volume. Since U(r,r') is localized in its spatial coordinates 
due to confinement of quarks and gluons, the potential receives finite volume effect 
only weakly in a large box. Therefore, once U is determined and is appropriately 
extrapolated to L — > oo, one may simply use the Schrodinger equation in the infinite 
space to calculate the scattering phase shifts and bound state spectra, which can be 
compared with experimental data. With this approach we successfully extract not 
only the nucleon potential Dl rt also the hyperon potential® below inelastic 

threshold in QCD. 

Although this method is shown to be quite successful in order to describe elastic 
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hadron interactions, the hadron interactions generally leads to inelastic scatterings 
as the total energy of the system increases. In this paper, we extend our method 
to such inelastic scatterings in order to extract hadron interactions in general. In 
Sec. [2j we briefly summarize our method previously used to extract the potential in 
the elastic scattering. In Sec. [3j we present our main idea to analyze the inelastic 
scattering from lattice QCD in the finite volume. We here discuss the scattering such 
as A + B^-C + D scattering, where A, B, C, D represent some 1-particle states. 
This is a simplified version of the baryon scattering in the strangeness S = — 2 and 
isospin 1 = channel, where AA, NS and ZJU appear as asymptotic states if the 
total energy is larger than 2m s- In Sec. HI we discuss the extension of our proposal 
to the scattering with particle productions such as A + B — > A + B + C . In Sec. 
we summarize this paper together with recent applications of our method in lattice 
QCD. A preliminary account of these results is given in RefJT]). 

§2. Hadron interactions below threshold: elastic scattering 

In this section, we summarize our strategy, previously used to extract the po- 
tential between two hadrons below inelastic threshold in lattice QCD. 

2.1. Nambu-Bethe-Salpeter wave function 

A key quantity in our method is the equal-time Nambu-Bethe-Salpter (NBS) 
amplitude, which we call the "NBS wave function" throughout this paper. Let us 
consider the following NBS wave function for a particle A and a particle B in QCD 
with total energy W in the center of mass system (i.e. the total three-momentum 
P = 0) in the infinite box; 

A(r)e- Wt = ]im(0\T{ip A (y,t + S)<p B (x,t)}\AB-W,P = 0) (2-1) 

where the relative coordinate is denoted as r = x — y. Here the local operators for 
the A and the B, which might be composite, are denoted by (fA(x,t) and ipB(y,t) 
with possible indices such as spinor or flavor being suppressed for simplicity. The 
QCD vacuum is denoted by |0), while the state \AB; W, P = 0) is a QCD eigenstate 
with the same quantum numbers as the AB system. Note that \AB; W, P = 0) can 
be taken as a product of 1-particle asymptotic state |-A)in <8> |-B)in in the infinite box, 
while this is not true in the finite box. 

If the total energy W = E A + = J m\ + k 2 + J m\ + k 2 is smaller than 

the inelastic threshold E t h, satisfies 

r ^3 jpr fiA , E A 

+ Z(r) (2-2) 

where pa = (E A ,p), Ua = {E^,k), ks = (Ej?,—k) are on-shell 4 momenta, while 
Pb = (W — E A , —p) is generally off-shell. Therefore T A b,AB (<Zi, 12, 13, Qi) is in gen- 
eral off-shell T-matrix, defined through the connected four-point Green's function 
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AB,AB (Pi iP2'iP3iPa) as 

g AB,Ab(pi>P2]P3,P4) = (2vr) 4 5 4 (pi +P2-P3-P4) 

x iD A (pi)iD B (p2)iT A B,AB ipi,P2',P3,PA) iD A (P3 ) (p$ • 3) 

where Da(b){p) 1S the free propagator for a particle -A(-B) in the momentum space, 
which does not contain the negative energy part (i.e. the contribution of the corre- 
sponding anti-particle). In the above expression, Z(r), which is exponentially sup- 
pressed for large r = \r\ as e~ cr with c oc E 2 h — W 2 > 0, represents contributions 
from other than elastic scattering A + B — > A + B. 
With the partial wave decomposition that 

4>J B (r) = 4tt ^AB,i(r, k)Y lm {n r )Y lm {Q h ) (2-4) 

l,m 

where k = \k\, Y\ m is the spherical harmonic function and fi r is the solid angle of 
the vector r, one can show for the large r that 

^ (r ,, hill !Ht^+«S, ( 2 . 5 ) 

where the "phase shift" 5i(k) is the phase of the S-matrix of the A + B — > A + B 
scattering for the partial wave I. Therefore the NBS wave function is indeed the 
"wave function" which describes the AB — > AB elastic scatteringP'l'^'® 

In the finite volume, restricted values of k denoted by k n can be realized to satisfy 
the boundary condition. From the energy of two particle Wl in the finite volume, 
one can determine the phase shift di(k n ) through Liischer's formula,^ where k n is 
determined from Wr = Et + E? . 

2.2. Strategy to define the potential 

In this subsection, we summarize our strategy to define the "potential" in QCD. 

(1) We choose the field operator tpA and (p B . If these operators are composite, 
there are many choices to create the same one-particle state. For example, we 
take the local operator for nucleon in our previous calculations for the nucleon 
potential^.!!® 

(2) We then measure the NBS wave function, defined by 

i>\ V B (r)= Urn (0\T{<p A {x + r, 5)<p B (x, 0)}\AB;W, P = 0). (2-6) 

(3) Motivated by the fact that the NBS wave function describes the elastic scatter- 
ing in the large r, we define the non-local potential as 

[E£ B - H$ B ] ^ b {x) = J d 3 y U(x;y)^ B (y), (2-7) 

k 2 V 2 

rpAB _ K ttAB _ _ v 

^AB Lp-AB 

where k = \k\ and the reduced mass [iab is defined by 1/hab = 1/mA + l/m B . 



4 



HAL QCD Collaboration 



(4) We then perform the velocity (or derivative) expansion that U(x; y) = V(x, V)5 3 (x- 
y). In the case of the NN scattering, for example, we have^ 

V{r, V) = V (r) + F ff (r)(«n ■ a 2 ) + V T ir)S l2 + V^{r)L ■ S + 0(V 2 ) (2-8) 

where <Xj is the spin operator of i-th. particle, S\ 2 is the tensor operator given 
by 

Si2 = ^i<Ti-r)i<T 2 -r)-i<T 1 - t T 2 ). (2-9) 

In the above expansion, the first three terms are of the leading order (LO), 
which do not contain derivatives, while the 4-th term is of the next leading 
order (NLO) with one derivative. 

(5) Once we obtain the potentials, we solve the Schrodinger equation with these 
potentials in the infinite volume to obtain physical observables such as the phase 
shift and binding energy. Note that the exact value of the scattering phase 
shift <5(/c) is obtained from this Schrodinger equation, while <5(/c') at k' ^ k 
is approximated one as long as the derivative expansion for the potential is 
truncated at the finite order. 

In order to extract the NBS wave function on the lattice, we evaluate the 4-point 
correlation function as 

g AB ix,y,t-t ;J p ) = (0\Ty A ix,t)ip B iy,t)}J AB it ;J P )\0) (2-10) 
= ^2A n (0\T{ip A ix,0) VB iy,0)}\W n )e- w ^-^ (2-11) 

n 

-)• A ipJ B ir = x-y; J p ) e - W ^\ t - t -)• 00 (2-12) 

where A n = {W n \J AB ito; J p )|0), and J AB ito; J p ) is some source operator which 
create 2-particle states of AB at to with fixed total angular momentum J and parity 

p[Ql,[2jl,[3l 

2.3. Frequently Asked Questions 

There exist some questions to the definition of the potential in the previous 
subsection. 

(1) Does the potential depend on the choice of operators ip A ,p B ? 

Yes. The potential of course depends on the operators (p A ,tp B from which the 
NBS wave function is defined. The choice of the operators can be regarded as 
the "scheme" to define the potential, since the potential itself is not a physical 
observable so that it can be scheme-dependent. If we calculate the physical 
observables, however, we obtain the unique result irrespective of the choice for 
the operators ip A ,ip B as long as Uix; y) is evaluated exactly. This is quite anal- 
ogous to the running coupling, which is of course scheme-dependent. Physical 
matrix elements do not depend on the scheme of the running coupling, as long 
as they are evaluated exactly. 

(2) Does the potential depend on the total energy at which the NBS wave function 
is defined ? 
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By construction, U(x;y) is non-local but energy- independent, while the defini- 
tion 

V k (x)=E k - ^ Jb{x) (2-13) 

gives a local but energy(momentum)-dependent potential. It is easy to see that 
Vfe(x) is equivalent to U(x,y) since the number of degrees of freedom of fc is 
equal to that of y in the finite volume. 31 From the fc-dependence of the Vk(x), 
we therefore can determine the higher order terms of the derivative expansion 
V(x, V). It turns out, however, that the fc-dependence of Vk{x) for the nucleon 
potential is very small in quenched QCD between k ~ MeV(£f N ~ MeV) 
and k ~ 240 MeV(£f N ~ 45 MeV) for our scheme, 13 ^, IS} .OS where the local 
composite field of three quarks is used for nucleon operator. 

§3. Hadron interactions above inelastic threshold 

3.1. NBS wave function in inelastic scattering 

We now consider the scatterings A + B — > A + B and A + B — > C + D. We 
assume that mA + tub < mc + < W, where W = E£ + is the total energy 

of the system with Eff = J m 2 x + fc 2 . In this situation, the QCD eigenstate with 
the quantum numbers of the AB state and center of mass energy W is expressed in 
general as 

\W) = c AB \AB, W) + c CD \CD, W) + --- (3-1) 
\AB,W) = \A,k) ia ®\B,-k) ia , \CD,W) = \C,q) ia ®\D,-q) ia , (3-2) 

where W = E£ + E% = E° + We define the following NBS wave functions, 

if>AB(r, fc) = lim {0\T{<p A (x + r, 8)<p B (x, 0)}\W), (3-3) 

<5— 5-0+ 

fcD(r,q)= lim {Q\T{yc{x + r,8)<p D (x, 0)}\W), (3-4) 

<5— 5-0+ 



ipAB(r,k) = \JZ A Z B 

(>) 3 p2 _ fc 2 _ ie mE A 



which can be expressed as 

[ ik . r f d 3 p Ep+Etrr.ABAB, , , , 

(3 



d 3 p e*- E^ + Ej D 
+ °CD J—^—^ T2 ~ AJjr „ A T {PA,PB,qc,qD) 



VcD{r,q) = yJZ c Z D 

d 3 p e ^ r E° + Ef 



d 3 p e ipr Ep + E° 



(3 
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where p A = (E A ,p), p B = (W - E$,-p), p c = (E°,p), p D = (W - E°,-p), 
k A = {E A , k),k B = {E B , -k), q c = {E%, q) and q D = (E° , -q). 
Introducing 

H ABAB { CD ){pMq)) = ^^T AB ^ CD \p A ,p B ,k A ,k B {qc,q D )) (3-7) 
H°W°)<p,k(q)) = ^^T c ^ AB ^\p c ,p D ,k A ,k B {q c ,q D )) (3-8) 



(3-9) 
(3-10) 
(3-11) 



and using the partial wave decomposition such that!*)! 

^XY(r,k) = 4ir ^ i l il> l XY (r, k)Y lm {Q r )Y lm {Q k ) 



l.m 



e ikr = ^J2i l Mkr)Y lm (f2 r )Y lm (f2 k ) 



l.m 



H XY > vz (p, k) = A^ H? Y ' VZ (p, k)Y lm (f2 p )Y lm (f2 k ) 



l.m 



with XY, VZ = AB or CD, we have 



^Ab(?~,&) = V z aZ b c ab lji(kr) + J Tj^^2 



+ CCD 



p 2 dp 



2tt 2 p 2 — k 2 — ie 



H. 



p z — k 2 — i< 

AB,CD 



H 



AB,AB 
I 



(p,k)ji(pr) 



{p,q)jl{pr) 



(3-12) 



^ l cD(r>q) = V z cZ D 

' p 2 dp 1 



CCD ^ji{qr) + J 



p 2 dp 1 
2tt 2 v 2 — q 2 



+ c AB 



-H 



CD.AB 



p z — g z — ie 
(P,k)ji(pr) 



-HF D,CD (p, q)ji(pr) 



j 2tt 2 p 2 — q 2 — ie 1 
Here the spherical harmonic function Y\ m is normalized as 



df2 r Yi m (f2 r )Yii m i(f2 r ) = 6u'8 r 



(3-13) 



(3-14) 



with the solid angle Q r of the vector r, and the spherical Bessel function ji(x) is 
given by 



ji(x) = (-x) 1 



1 d \ l ( sinx\ sin(x — Itt/2) 



x dx 



x — > CO. 



(3-15) 



The p integral givesP®' 1 ^ 



^AB( r 5 fc ) = V Z aZb 



'\\b <i jl(kr) + -^Hf B ' AB (k, k) {ni(kr) + iji(kr)} 



*' Here we ignore spins for simplicity. 
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+ ccd ^Hf B ' CD (k, g) { ni (kr) + ij, (At)} 



+ l l AB (r) (3-16) 
^odO", Q) = V z cZ D ccd \ji(qr) + -^-H^ D ' CD (q, q) {n^qr) + ij t (qr)}} 



(3-17) 



where X^-y(r), which represents all contributions except those from the pole at p = k 
or p = q, is exponentially suppressed for large r, and another spherical Bessel function 
n/ (x) is given by 



m{x) = (-x) 1 

The unitarity relation 



l_d_ 

x dx 



^ cos x ^ 



cos(x — Itt/2) 



x — ¥ oo. 



X 



T - T f = 7T f r 

for the on-shell T-matrix T of the 2 channel scattering gives 



(3-18) 



(3-19) 



Pi 



0{W) 



e i5 K w ) sin 5l(W) 

e iS ?W sin 6f(W) 



cos 9(W) -sin6(W) 
sin0(W) cos 9(W) 



7=1,2 J = l,2 



O" 1 (VF)(3-20) 
(3-21) 



where tfj(W) is the scattering phase shift , whereas 9{W) is the mixing angle between 
1 and 2. Here index 1 represents AB while 2 represents CD. With this notation 
pi = k and P2 = q- We therefore obtain 



Pi TjI,J( \ 

^tf,' (P/,PJ) 



O(W) 



e iS ?^ sin 5f(W) 



0~\W) 



I, J 
3-22) 



Using these results, the NBS wave functions of the 2 channel system behave for 
large r as 



ip l AB {r,k) 
%d{ 



jl(kr) 

ji{qr) 



CAB 
CCD 



+ 



ni(kr) + iji(kr) 

ni(qr) + iji(qr) 



o<w)( ^ iSHw) ^Uw) o 

'> e^Wsin^(^) 



Q~ 1 {W) 



CAB 
CCD 



3-23) 



where W X y = V'xy/V-^x-^y- This expression shows that the NBS wave function 
for large r agree with scattering waves described by two scattering phases 8\{W) 
(i = 1, 2) and one mixing angle 0(W). 



8 



HAL QCD Collaboration 



3.2. Coupled channel potentials 

Let us now consider QCD in the finite volume V. In the finite volume, \AB, W) 
and \CD,W) are no longer eigenstates of the hamiltonian. True eigenvalues are 
shifted from W to Wi = W +0(V~ l ) (i = 1, 2). By diagonalization method in lattice 
QCD simulations, it is relatively easy to determine W\ and W%. With these values 
Liischer's finite volume formula gives two conditions, which, however, are insufficient 
to determine three observables, 5f, 5f and 6. fSee fT7|l . [T8]) . [T9|) for recent proposals 
to overcome this difficulty.) We here propose alternative approach to extract three 
observables, 5j, 5f and 6, in lattice QCD through the above NBS wave functions. 
We consider the (normalized) NBS wave functions at two different values of energy, 
W\ and W 2 , in the finite volume: 

VUfl(r,fci) = J— lim .(0\T{<p A (x + r,S)<p B (x,0)}\Wi) (3-24) 
\J Z A Z B <5^0+ 

V>cr>(r,g f )= -=^= lim (0\T{<p c (x + r,5)<p D (x,0)} \ Wi), i = 1,2. (3-25) 

(We here omit ~ on tp.) We then define the coupled channel non-local potentials from 
the coupled channel Schrodinger equation as 

[E^ B -H^ B ] ip A B(x,ki) = J d 3 y U AB ,AB(x;y) ^ AB (y,ki) + J d 3 y Uab,cd{x;v) ipCD(y,Qi) 

(3-26) 

[Eg D - Hg D ] i) C D{x,ki) = j d 3 y U CD ,AB(x;y) ip AB (y,ki) + J d 3 y U CD ,CD(x;y) ^CD(y,qi) 

(3-27) 

for i = 1,2. As before we introduce the derivative expansion as 

U X Y,vz(x;y) = V XY ,vz{x,V)5 3 (x - y) = [V XY ,vz(x) + O(V)] 5 3 (x - yp-28) 
and at the leading order of the expansion, we have 

K AB (x,ki) = [e£ b - Hq B ] ip AB {x,ki) = V ABi ab{x) if> AB (x,ki) + V ABi cd(x) ^00(33,^) 

(3-29) 

K C D{x,qi) = [E^ D - H$ D ] ipCD{x,ki) = V C d,ab(x) ip AB (x,ki) + U C d,cd(x) ip C D{x,qi). 



(3-30) 



These equations for i = 1, 2 can be solved as 

Vab,ab(x) Vab,cd{x) \ _ ( K AB (x,k{) K AB (x,k 2 ) 

V C d,ab(x) Vcd,cd{x) ) \ K G D{x,qi) K C D(x,q 2 ) 

tl> AB (x,ki) ip AB (x,k 2 ) x 

^Cd(«, <Zi) ipcD(x,q 2 ) 



(3-31) 



Once we obtain the coupled channel local potentials Vxy,vz(x), we solve the cou- 
pled channel Scrodinger equation in infinite volume with some appropriate boundary 
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condition such that the incoming wave has a definite I and consists of the AB state 
only , in order to extract three observables for each /, 5j(W), 5f (W) and 0(W), at all 
values of W . Of course, since Vxy,vz is the leading order approximation in the veloc- 
ity expansion of Uxy,vz{x; y), results for three observables 5]{W), 6f(W) and 9(W) 
atW^ Wi, W2 are also approximate ones and might be different from the exact val- 
ues. By performing an additional extraction of Vxy,vz(x) at (W3, W4) 7^ (Wi, W2), 
we can test how good the leading order approximation is, as in the case of the elastic 
scattering.^ 

The above result can be easily extended to the coupled channel among n states 
AjBj (I = 1,2, ■•■ , n) . Local potentials at the LO in the velocity expansion are 
given by 

Vi,j( r ) = E Kjir^i^X-'i^Mj (3-32) 

M=l 

Kj(r, *4) = [E r M - Hi] Mr, *4) (3-33) 
where X~ l (r) is the inverse of the n x n NBS wave function matrix 

X(r)jj = V/(r, fc$) = —L== lim (0\T{tp Al (x + r,8)tp Bl (x,0)}\Wj)(3-34) 



the momentum kj satisfies the relation that Wj = y (kj) 2 + vn Aj + y (kj) 2 + m? Bl 
for J = 1,2, • • • ,n, Ej = {k 1 ^ 2 / {2^ AlBl ) and H* Q = -V 2 /(2// AjBj ). 

3.3. Inelasticity and non-locality 

We now consider the 2-channel problem again. Let us assume that the leading 
order approximation for the coupled channel potential works reasonably well. In- 
stead of considering the coupled channel potential, the effective potential for the AB 
channel is given by 

Uf ByAB (x,y) = V ABAB {x)5^{x - y) 

+ V ab , C d(x) e c D _ R cD _ cd 0> v)Vcd,ab(v), (3-35) 

where the non-locality becomes manifest in the second term. The magnitude of 
momentum q for the CD channel is expressed as 

4W 2 q 2 = {W 2 - (m c + m D ) 2 )(W 2 - (m C - m D f). (3-36) 

To estimate the magnitude of non- locality in eq. (|3-35p . we here ignore the Vcd, CD 
term in the denominator. In this case, if the total energy W is below the inelastic 
threshold mc + mo such that q 2 = —M 2 < 0, we have 

1 (x,y) = - 2flCD - e- M \ x -V\ (3-37) 



Hq D — E^ D ' y; 4w\x-y, 

Therefore, non-locality of the potential U A S B AB (x,y) is exponentially suppressed as 
long as M is large enough. As W approaches mc -\-mjj, however, non-locality of the 
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potential becomes larger and manifest. In this simple example, this non-locality of 
the single channel potential can be completely removed by introducing the CD state 
and coupled channel potentials between AB and CD. Therefore, in practice, the 
nondocality caused by inelastic final states whose thresholds are closed to the total 
energy W is expected to become milder for the coupled channel potentials including 
these states. 



§4. Extension: Inelastic scattering with particle production 

The method considered in the previous section can be generalized to inelastic 
scattering where a number of particles is not conserved. For illustration, we consider 
a case that the scatterings A + B — > A + B and A + B — > A + B + C occur and the 
total energy W satisfies tua + m B + mc < W < tua + ms + 2mc- 

We consider the following NBS wave functions in the center of mass system: 

1>ab{*) = -^=^(0\^A(r + x,5)ip B (r,0)\W) (4-1) 

V Za^B 5^0+ 

^ BC (x, y) = 1 lim ,(0\<p A (r + x + 8)<p B (r + y, 0)^(r, -<5)|(t¥£) 



where 



\W) = Ci |fc) in <8> | - fc) in + C 2 \q x )hi <8> \q y ~ q x ^^)m ® \ ~ Qy ~ Qx^-}in(4-3) 



with 



W = k 2 + m\ + y k 2 + m| 

= y/<ll + < + sj^y ~ ^? + m% + yj{q y + q x ^Y^i (44) 

and 1/ fiBC = 1/mB + 1/mc- Here y = r B — rc is a relative coordinate between B 
and C with the reduced mass fi B c, while x = va — ~R.bc is the one between A and 
the center of mass of B and C with R B c = {wibTb + rncrc)/(rnB +mc)- We here 
assume the property that i^abc nas asymptotic behavior of the scattering wave of 
A + B + C as \x\, \y | — > oo, as shown for ip^B * n ^ ne previous section. Although this 
property has not been shown so far, it is reasonable to assume this. We leave the 
proof for this important property to the future investigation. 

We define the nondocal potential from the coupled channel equations as 

Kf B (x) [E£ B - H$] ^ B {x) = Jd 3 z U A b,ab(x-z)^ b {z) 

+ J d 3 z d 3 w U A b,abc(x;z,w)iPab C (z,w) (< 
Kf BC (x, y) [Efj BC + Ef y c - H*f c - H™] ^ BC (x, y) = J d 3 z U A bc,ab(x, y 
x ^ab{ z ) + J d 3 z d 3 w U A BC,ABc(x,y;z,w)^ BC (z,w) 



Hadron Interactions above Inelastic Threshold 



11 



ttAB 
n 0,x — 




u a,bc 

n 0,x — 




ttBC 

H 0,y ~ 


-V 2 
v y 




2/J>A,BC ' 






k 2 


E A,BC = 




P BC 




I^ab' 


2[J>A,BC ' 


IHBC 



where 

(4-6) 
(4-7) 

with another reduced mass defined by 1/ ' ha,bc = + V( m £ + m c)- 

We consider the following velocity expansions 

Uab,ab{x- z ) = [Vab,ab(x) + 0(V X )} 5 3 {x - z) (4-8) 
U A b,abc(x; z, w) = [V AB ,abc(x, w) + 0(V X )] <5 3 (x - z) (4-9) 
U A bc,ab(x, y, z) = [V A bc,ab(x, y) + 0(V X )] 5 3 (x - z) (4-10) 
U abc, abc { x ! y; 2, to) = [Vabc,abc(«, 2/) + 0(V X , V„)] -5 3 (a; - z)5 3 (y - u#ll) 

where the hermiticity of the non-local potentials gives Vab,abc( x , y) = Vabc,ab(x, y)- 

At the leading order of the velocity expansions, the coupled channel equations 
become 

K%j(x) = V A b,ab( x )^b( x ) + Jd 3 w V A b,abc(x,w)^ bc (x,w) (4-12) 

K A bc(x, y) = V A bc,ab(x, y^^ B (x) + V A bc,abc( x , y)^)% BC (x, y). (4-13) 

By considering two values of energy such that W = W±,W2, we can determine 
Vabc ab and Vabc,abc from the second equation as 



( V AB c,AB(x,y) V AB c,ABc(x,y) ) = ( K™ BC (x,y) K^ c (x,y) ) 

^ B ( X ) ^ B (x) y\ 

^ BC (x,y) ^ BC (x,y) J { ) 

Using the hermiticity Vab,abc{ x ,v) = VABC,AB( x ,y), we can extract Vab,ab from 
the first equation as 



Vab,ab(x) 



n B (x) L 



k ab( x ) - J d 3 w V A BC,AB( x ,w)tpJ BC (x,w] 



(4-15) 



for W = Wi,W2- A difference of Vab,ab( x ) between two estimates at W\ and W2 
gives an estimate for higher order contributions in the velocity expansions. 

Once we obtain Vab,ab, Vab,abc = Vabc,ab and Vabc, abc-, we can solve the 
coupled channel Schrodinger equations in the infinite volume, in order to extract 
physical observables. As W increases and becomes larger than + m B + nmc, 
the inelastic scattering A + B—tA + B + nC becomes possible. As in the case of 
A + B — > A + B + C in the above, we can define the coupled channel potentials 
including this channel, though calculations of the NBS wave functions for multi- 
hadron operators become more and more difficult in practice. 
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§5. Summary 

In this paper, we propose extensions of the method of extracting potentials 
through NBS wave functions to the case where inelastic scatterings becomes im- 
portant. We first consider the case that A + B — > C + D scattering occurs and 
present an explicit formula to extract the coupled channel potentials between AB 
and CD. The general formula for the scattering among n states is also given. An 
extension to the case where the particle production occurs during the scattering such 
as A-\- B — > A + B + C is also considered. This can also be extended to more general 
cases such asA + B^-A + B + nC (n = 1, 2, 3, • • • .). 

Recently the potential method has been used to extract potentials for the flavor 
SU(3) limit in lattice QCD where up, down and strange quark masses are all equalP^ 
In this limit there exist 6-independent potentials corresponding to the irreducible 
representations of the flavor SU(3), and among these, the flavor singlet potential is 
strong attractive, suggesting an existence of a bound state in this channel. This 
expectation has been confirmed by more detailed study in lattice QCD and a bound 
state of two up quarks, two down quarks and two strange quarks, called the H 
dibaryon, indeed exists in the flavor SU(3) limit!^ To study the property of the 
H dibaryon in the real world where the strange quark is much heavier than up and 
down quarks so that the flavor SU(3) symmetry is broken, the extension presented 
in this paper for the coupled channels is indeed necessary: The SU(3) breaking 
in nature appears in the octet baryon mass as mjv = 939 MeV, m\ = 1116MeV, 
ms = 1193 MeV and m= = 1318 MeV. Therefore thresholds of two baryon systems 
with strangeness S = —2 and ispspin / = 0, to which the H dibaryon belongs, are 
given by 

W A a = 2232 MeV < W NE = 2257 MeV < W ss = 2386 MeV. (5-1) 

A study for coupled channel potentials in this case has started using the 2+1 flavor 
lattice QCD and a preliminary result has already been presented.^ More details 
results for this study will be published soon. 
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